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Monitoring and visualizing specimens at a large penetration depth is a challenge. At depths of hundreds of 
microns, several physical effects (such as, scattering, PSF distortion and noise) deteriorate the image quality 
and prohibit a detailed study of key biological phenomena. In this study, we use a Bessel-like beam 
in-conjugation with an orthogonal detection system to achieve depth imaging. A Bessel-like penetrating 
diffractionless beam is generated by engineering the back- aperture of the excitation objective. The proposed 
excitation scheme allows continuous scanning by simply translating the detection PSF. This type of imaging 
system is beneficial for obtaining depth information from any desired specimen layer, including 
nano-particle tracking in thick tissue. As demonstrated by imaging the fluorescent polymer- tagged- CaC0 3 
particles and yeast cells in a tissue-like gel-matrix, the system offers a penetration depth that extends up to 
650 jum. This achievement will advance the field of fluorescence imaging and deep nano-particle tracking. 

Fluorescence imaging techniques have dramatically changed biology and medicine. Recently, it was proved 
beyond a doubt that, fluorescence imaging techniques are capable of going beyond the classical resolution 
limit 1 , thereby revealing molecular dynamics with nanoscale resolution. Broadly, these super-resolution 
techniques may be categorized based on either point spread function (PSF) engineering or photophysical prop- 
erties of the molecule. Techniques using PSF engineering include, 4PI 2,3 , stimulated emission depletion (STED) 
microscopy 4,5 , aperture engineering/multiple excitation spot optical (MESO) microscopy 6,7 , spatially structured 
light illumination 8 , non-linear patterned excitation microscopy 9 , and standing-wave total internal reflection 
fluorescence (TIRF) microscopy 10 in addition to others 1112 . Photophysical property-based super- resolution tech- 
niques include photo -activation localization microscopy (PALM) 13 , stochastic optical reconstruction microscopy 
(STORM) 14 , fPALM 15 , ground-state depletion (GSDIM) microscopy 16 and IML-SPIM 17 . Other techniques that do 
not belong to these categories are evanescent field-based techniques such as, TIRF 1819 and HILO 20 . Notably, 3D 
image reconstruction in conjunction with these super- resolution techniques produce images of the highest 
quality 5,21 " 23 . PSF engineering utilizes optical techniques, such as, interference and spatial filtering to shrink 
the system PSF, whereas localization techniques use photophysical phenomena, such as, photoactivation and 
depletion, to achieve super-resolution. The down side of these photophysical techniques is the sacrifice of 
temporal resolution to gain spatial resolution because these techniques collect thousands of frames to construct 
a single super-resolution image. Nevertheless, the contribution of some of the super- resolution techniques are 
noteworthy and considered to be a tool for future development. The contributions include, resolving lateral 
elements of the synaptonemal complex and the detection of chromosomes as twisted strings 24 . Localization-based 
microscopy has revealed the clustering nature of chemotaxis receptors in the E. Coli cell membrane, which 
supports the fact that stochastic self-assembly is the cause behind periodic distribution of the receptors 25 . The 
study of vesicle dynamics reveals the release of lipid probes (from granules) into the membrane without fusing 
thereby revealing a new mechanism 26 ' 27 . Multicolor TIRF microscopy is extremely useful for observing cellular 
dynamics near the cell membrane 28 . Diffraction-unlimited optical data storage has been achieved with poly- 
chromic GFP using STED imaging 29 . A variant of STED i.e., iso-STED, was able to resolve and study the 
distribution of mitochondrial proteins 30 . It was observed that, nuclear pore complexes are adjoined in peripheral 
and lamin heterochromatin using structured illumination microscopy 31 . There are many more studies revealing 
the benefits of super-resolution microscopy ranging from biology to optics. Unfortunately, the super-resolution 
capabilities are limited to a depth of < 150 microns, which is largely attributed to scattering in complex biological 
specimens. 
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Figure 1 | Schematic diagram of the proposed imaging system, in which a Bessel-like beam illuminates the sample with theta-detection at 90° . Spatial 
filtering using a binary optical mask behind the illumination objective results in the generation of a Bessel-like beam. Position 1 and position 2 of the 
detection arm represents two different scanning positions along the optical z- axis for depth imaging. 



Most of these techniques have impressive super-resolution cap- 
abilities but are limited in terms of depth imaging and suffer from 
poor axial resolution. Another limitation is the point-by-point-based 
slow scanning employed in most of the existing fluorescence imaging 
systems. A Bessel beam-based illumination scheme has shown prom- 
ise for monitoring and imaging large specimens 32 . Because a pencil- 
like region is illuminated, the specimen is well-protected from photo - 
bleaching in comparison to confocal microscopy which illuminates 
the layers above and below the optical plane. We plan to lift some of 
these restrictions using spatial filtering techniques by employing a 
Bessel-like diffraction-less beam for excitation in- conjugation with 
theta-based orthogonal detection scheme 33 ' 34 . A Bessel-like excitation 
beam is produced by placing a binary mask at the back- aperture of 
the excitation objective. The generation of a Bessel-like beam using 
an aperture engineering technique has been reported in the litera- 
ture 35 " 37 . Notably, the proposed aperture engineering technique sub- 
stantially reduces the ringing effect, thereby enabling fine sampling 
of the specimen. Reducing the ringing effect in a Bessel-like beam mi- 
nimizes the background fluorescence, thereby improving the signal- 
to-background ratio. Such an imaging technique allows continuous 
scanning of specimens along the z-axis at variable depths. This tech- 
nique substantially simplifies the imaging system and eliminates the 
need for a specialized illumination technique such as structured 
illumination. The advantage of the proposed technique is its simpli- 
city and adaptability to existing imaging techniques. Specifically, this 
technique is well-suited for tracking functional nano -particles at 
greater penetration depths, enabling better understanding of com- 
plex biological processes deep inside large specimens such as, tissue, 
membrane and even live embryos. 

Imaging deep inside thick biological tissue is severely hampered by 
photon scattering, which results in intensity modulations and phase 
mis-matches. Most fluorescence imaging techniques use a Gaussian 
beam for illumination. The Gaussian illumination technique is well- 
established and useful for shallow imaging (less than 150 microns), 
but fails to produce quality images at larger depths (few hundreds of 
microns). Two-photon excitation has been specifically introduced to 
overcome the above drawbacks. However, the output fluorescence 
is severely compromised as a result of the small absorption cross- 
section, thereby limiting its use for fluorescence imaging 38 . This calls 
for a beam that maintains its shape and size at larger penetration 
depths. Such a diffraction-less beam is known as a Bessel beam. It is 
the self-reconstruction property of the Bessel beam that maintains 
the profile of the beam in thick scattering specimens 39 . 



Results 

In this paper, we demonstrate the generation of a Bessel-like beam 
and show that it can reach deep inside the sample. Another aspect of 
the generated Bessel-like beam is the reduction of ringing effects that 
enhance background fluorescence. By using an orthogonal detection 
system, we demonstrate the scanning capability of the imaging sys- 
tem with a considerable reduction in the background fluorescence. 
This study shows an elongated Bessel-like beam and larger penetra- 
tion depths. 

System PSF and its characterization. Studies were carried out to 
characterize the PSF of the imaging system. The excitation and emis- 
sion wavelengths were chosen as 532 nm and 560 nm respectively. 
The schematic diagram of the proposed imaging system is shown in 
figure 1. The imaging system consists of two arms: the excitation arm 
is basically an aperture-engineered objective lens. The spatial filter 
belongs to the family of binary filter that allows light from the 
peripheral regions of the spatial filter. Spatial filtering followed by 
focusing using an objective lens results in a Bessel-like penetrating 
beam whose dimension (along the optical axis) can be tuned (by the 
stop angle 0 Y or equivalently by the radius r of the spatial filter) to 
generate a uniform field strength throughout the depth-of-focus. 
This technique simultaneously excites a pencil-like region encom- 
passing all of the specimen layers and the detection system is placed 
orthogonal to the optical z-axis, thereby enabling quick scanning 
and high axial resolution. In this configuration, the axial resolution 
is essentially the lateral resolution of detection objective. The de- 
tection objective (0 2 ) can be rapidly scanned (along the optical z- 
axis) to obtain structural information throughout the specimen 
layers (see figure 1). 

The system PSF at an excitation (oc exc ) and detection (a det ) semi- 
aperture angle of 60° is shown in figure 2. The lateral XY- plane shows 
the ringing effects associated with Bessel beams. This effect disap- 
pears in the system PSF (X7-plane), thereby producing an ellipse-like 
lateral PSF. The z-profile of the Bessel-like beam (see, YZ and XZ 
planes) is shown in figure 2. Excitation PSF multiplied by the detec- 
tion PSF results in a high resolution spot along the optical z-axis. The 
FWHM of the PSF along the z-axis is approximately 300 nm below 
the classical diffraction limit. This is impressive considering poor 
axial resolution of the existing fluorescence imaging system. 

Depth scanning. Fast scanning at a depth of hundreds of microns is 
essential for obtaining vital information about a specimen. In this 
paper, we propose Bessel beam-based depth imaging to gain depth 
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Figure 2 | The excitation, detection and system PSF at an excitation and detection aperture angle of a m = 60° and a det =60°, respectively. Linearly 
polarized light of wavelength 532 nm is employed for excitation and the fluorescence is observed at a wavelength of 560 nm. The resulting axial and lateral 
beam profile for the generated Bessel-like beam (central lobe) along with the side rings is clearly visible. The system PSF are obtained by multiplying the 
excitation and detection PSFs. Scale bar = 1 \im. 



information. Figure 3 shows the system PSF at various stages of the z- 
scanning process starting from z = 0 to z = 1.92 fim. Specifically, the 
system PSFs are shown for Az = 0 nm, 510 nm, and 1.47 fim. The 
PSF retains its shape and size throughout the scanning process along 



the optical axis, but this may not be true with a real sample which 
is non-isotropic and causes prominent scattering. Nevertheless, this 
shape retention is beneficial for depth imaging applications. Because 
of the non-diffracting nature of Bessel-like beams, the technique may 
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Figure 3 | The PSFs obtained by performing z-scans along the optical axis. The detection objective is scanned along the optical axis to obtain data at 
varying sample depths: Az = 0, Az = 510 nm and Az = 1470 nm. In this figure, Az represents the shift from the origin as depicted schematically in 
Figure 1. At varying depths, the system PSF preserves its shape and size. The effects of prominent side rings are clearly minimal. Scale bar = 1 /urn. 
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Figure 4 | PSF resizing by varying the transmission angle AO — a — 0 x -> where the stop angle & x is determined by the spatial filter. The excitation PSF 
axially elongates with decrease in AO and ranges from 1.56 fim to > 3.84 /urn. For large AO, the excitation PSF resembles the PSF obtained using plane wave 
illumination. The intensities are normalized for each image. Scale bar = 1 /urn. 



also perform well in scattering samples. However, emitted photons 
may undergo significant scattering as they travel through the spe- 
cimen along the orthogonal direction to reach the detector. Overall, 
the proposed system benefits as far as high resolution depth imaging 
is concerned and provides the flexibility to scan the desired specimen 
layer. Next, we show the effect of the filter characteristics on the 
system PSF. Figure 4 shows the excitation PSF for varying filter 
stop angles. The PSF scales with the stop angle (0i), as the trans- 
mission window AO = (ol — 0i) is varied from 5° to 20°. The cor- 
responding FWHM along the optical axis are FWHM Z = 1.56 /urn 
(for AO = 20°), FWHM Z = 2.94 \im (for AO = 10°) and FWHM Z > 
3.84 fim (for AO = 5°). Additionally, the FWHM Z becomes hundreds 
of microns for AO < 5°. Ideally, the system PSF resembles widefield 
PSF for small stop angles (i.e., large AO), resembling full aperture 
uniform illumination. This scaling makes the system flexible for 
many biological applications that demand PSF of specific dimen- 
sions in real time, especially when imaging large biological samples 
such as an embryo. Finally, the polarization- dependence of the 
excitation PSF is analyzed. Figure 5 shows, respectively the contour 
plots of the excitation PSF for both x-polarized and randomly po- 
larized light illumination. For randomly polarized light, a uniform 
decrease is observed in the field strength as the distance form the 
geometrical focus i.e, (x,y, z) = (0, 0, 0) increases along the transverse 
plane. This uniform decrease is completely different from that of 
linearly polarized light illumination, which shows elliptical field 
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structures as the distance from the focus increases along the y-axis. 
The field strength does not decrease significantly along the )/-axis, 
even as far as 1.92 \im from the geometrical focus. 

Experimental results. Following the computational study, we con- 
structed an optical system to generate a Bessel-like beam for exciting 
the fluorescent gel sample. The self- reconstruction property of Bessel- 
like beams facilitate navigations through thick and turbid media. The 
interference of the diffracted beam emerging from the engineered 
aperture mask results in the formation of Bessel-like beam. Aperture 
engineering has been quite successful for producing Bessel beams of 
desired dimension. Altering the filter parameters (aperture width 2r or 
equivalently, stop angle 6{) of the spatial filter mask varies the size in a 
simple way 37 . This adaptability is the advantage over axicon-based 
Bessel beam generation, where there is no tunability in the axial di- 
mension of the Bessel beam. 

The schematic diagram for obtaining the Bessel beam-based high 
resolution imaging is shown in figure 1. The input Gaussian beam is 
appropriately expanded to fill the back- aperture of the objective lens 
and subsequently subjected to the binary spatial filter before entering 
the illumination objective lens 0^ The diffracted waves undergo 
constructive interference along the optical axis, thereby forming a 
penetrative diffraction-less Bessel-like beam. The XY, XZ and YZ 
planes of the Bessel beam are shown in figure 6. The measurements 
were taken by translating the camera that was placed directly in the 
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Figure 5 | Polarization effect on the excitation PSF: (a) Linearly-polarized light illumination (polarized along x-axis), and (b) Randomly polarized light 
illumination. The contour plots clearly show an isotropic field distribution for randomly polarized light, whereas the field has anisotropy along the y-axis. 
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Figure 6 | Experimentally obtained illumination PSF along the XY, XZ and YZ planes showing the extended depth-of- field. A camera was placed 
directly in the beam path near the focal plane of the illumination objective and several images of its transverse profile were captured by repeatedly 
translating the camera through a known distance to obtain a 3D image. These orthogonal views are obtained from the 3D data. The Bessel beam spatially 
extends approximately 750 /urn along the optical axis. Intensity plots are shown along the x y y and z axes for detailed visualization of the field structure at 
and near the focal plane. The field is uniformly stretched along the optical z-axis showing the extended depth of the Bessel beam. 



beam path. The Bessel beam extends over a length of approximately 
750 [im (FWHM Z ) with a lateral (y) dimension FWHM y of approxi- 
mately, 50 /urn. These dimensions proves that, the generated Bessel- 
like beam is uniform over a fairly large axial distance, which is bene- 
ficial for precise microscopy imaging and spectroscopy deep inside 
the specimen. Moreover, the extent of the reported beam is larger 
than that earlier reported in the literature 32 . The advantages 
of using the proposed spatial filter instead of axicons are the easy 
tunability of the Bessel beam, the depth imaging capability and the 
low-cost 37 . 

The detection arm consists of a separate objective for collecting the 
fluorescence emanating from the fluorescent gel-matrix sample (see, 
sample preparation section for the preparation details of the sample). 
The detection arm is placed at an orthogonal angle offering high 
numerical aperture (NA) at long working distance apart from achiev- 
ing maximum elimination of incident light. After emerging from the 
back of objective lens 0 2 , the parallel rays are then filtered using a 
long-pass filter (with a cut-off wavelength of 550 nm) to remove the 
stray light and are finally focused on to the CCD camera. It may be 
reminded that, the system PSF is the product of excitation PSF and 
the detection PSF. In this system, the excitation PSF is the diffrac- 
tionless Bessel-like beam penetrating through all of the specimen 
layers, and the detection PSF is a dot-like structure (formed by the 
objective lens 0 2 ) which is placed orthogonal to the excitation object- 
ive (see, figure 1 and figure 2). The schematic diagram (see, figure 7A) 
shows the translation of detection PSF along the Bessel beam, thereby 
scanning the specimen along the length of the excitation PSF. The 
corresponding contour plots of the system PSF at varying depths of 
the fluorescent gel matrix are shown in figure 7B. The system PSF is 
obtained at a depth ranging from 50 fim to almost 1 mm. This range 
demonstrates the incredible penetrating property of the excitation 
PSF. Moreover, the system PSF preserves its shape to a depth of 650 
[im. The axial length (measured as the FWHM along the optical axis) 
and lateral (along transverse y-axis) extent of the Bessel-like beam are 
approximately 300 \im and 25 \im, respectively. The length of the 
Bessel beam can be varied by changing the diameter 2r of the spatial 
filter which is equivalent to changing the stop angle (0i). The adapt- 
able length of Bessel-like beam which depends on the spatial filter 
radius (shown in figure 1) provides the tunability required for vari- 
ous applications in bioimaging. 



We imaged coated CaC0 3 particle clusters with dimensions of 3 - 
10 fim in a gel matrix. In the present experimental setup, we have 
used both 10X, NA = 0.3 and 20X, NA = 0.4 objectives. Specifically, 
the imaging of CaC0 3 particles in the gel-matrix were carried out 
using 2 OX detection objective. The images were acquired at varying 
depths, from which only 2 images at a depths of 108 \im and 358 \im 
(shown by red circles) are shown in figure 8. The zoomed image of 
coated CaC0 3 -cluster is also shown for visualizing the structures. The 
image quality is impressive at a depth of 358 \im considering the fact 
that the gel matrix is highly scattering. The internal structure of the 
CaC0 3 cluster is clearly visible even at a depth of 350 \im and thus 
demonstrates in-depth high resolution imaging in a scattering med- 
ium (gel-matrix). This enhanced depth is a result of the self- recon- 
struction property of the Bessel beam 39 . Traces of the larger cluster 
are also visible in the image at position 1 (green circle) which are 
caused by the stray light and the large field- of- view of the detection 
objective. 

Imaging cells deep inside a tissue-like scattering specimen is of 
potential interest for understanding complex biological phenom- 
enon. In the present study, yeast cells were imaged in a tissue-like 
scattering gel-matrix. This type of specimen is important for asses- 
sing imaging capabilities in a scattering specimen while isolating 
other effects such as, auto-fluorescence and photobleaching. The 
yeast cells were coated with three-layers of fluorescently- tagged poly- 
mers. The coating process and sample preparation are described in 
Supplementary 1. Figure 9 shows the cell images at varying depths, 
Az = 296 fim (A), Az = 478 \im (B) and Az = 616 fim (C) (shown by 
the white arrow in figure 9). Similar to CaC0 3 imaging (in figure 8), 
traces of cells from other layers are also visible. Again, these traces are 
caused by stray light and the large field- of- view of the detection 
objective (0 2 ). Practically, the working distance for large NA object- 
ive (NA > 0.4), is usually very small (approximately few millimeters). 
Therefore, the optical alignment and data collection in the ortho- 
gonal detection configuration become very difficult. Moreover, the 
full-width half-maxima (FWHM) is calculated for several yeast sam- 
ples at varying depths as shown in figure 9. The intensity plots are 
shown in figure 9 (D,E,F) and the corresponding FWHM is found to 
be 2.04 fim, 2.36 fim and 2.12 \im at depths of Az = 296 \im, 478 fim 
and 616 \im respectively. These findings reveal that, the the FWHM is 
approximately preserved throughout the sample depth of 616 /im, 
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Figure 7 | (A) Cartoon showing the system PSF as it enters the fluorescent gel sample at varying z-depths along with the theta-detection geometry, 
(B) The experimentally obtained contour plots of the PSF at varying depth (ranging from 50 jum to 1.25 mm) of the fluorescent gel sample. The depth 
information is indicated on each image and denotes the distance of the right-most corner of each image from the coverslip. The depth penetration 
capability of the Bessel beam is quite evident. 



thereby demonstrating the diffractionless property of the Bessel 
beam. 

Finally, we performed depth-imaging studies with a state-of-art 
laser scanning confocal microscope (Zeiss, LSM 700) for comparison 
with the proposed technique. We used a confocal pinhole of 46.5 fim 
for obtaining the sectional images at varying specimen depths. The 
specimen comprised of fluorescently- coated Yeast cells encaged in 
a gel-matrix. Figure 10(a,b,c) shows the XY-, YZ- and ZX- sections of 
the specimen along with a cartoon (figure lOd) describing the 
imaging parameters (such as, the position of the coverslip, and the 
coordinates of the sample plane). From the YZ-plane (figure 10b), it 
is clear that the penetration depth of the confocal system is approxi- 
mately 227 /urn. A penetration depth of about 200 \im for a confocal 
system was recently reported in the literature 40 . Figures lOe and lOf 
show the magnified image (obtained using a PLANNeofluor 20X/0.4 
NA objective lens) of the specimen at a depth of 16 fim for direct 
comparison with the images at a depth of 227 fim (see, inset of 
figure 10a). To further demonstrate the system resolution, FWHM 
was obtained at varying depths. At depths of 16 fim and 227 ^m, the 
FWHM is found to be 2.29 /urn and 3.15 \im respectively (see, 
figure lOg and lOh). The broadening of the FWHM is clearly evid- 
ent, revealing the deteriorating resolution of the confocal system at 



greater depths. In contrast, the proposed Bessel beam-based imaging 
system retains the FWHM of the system PSF at much greater depths 
(see, figure 9). This gives an approximate estimate of the resolution of 
the proposed imaging system and its comparison with the confocal 
system. Similar images were also obtained for another specimen 
(fluorescent polymer- coated CaC0 3 particles embedded in the gel- 
matrix) as shown in the supplementary material (see, supplementary 
1, figure 4). 

Discussions 

The aperture engineering technique offers an efficient way of gen- 
erating a Bessel-like beam and the ability to reach a depth of hun- 
dreds of microns. The diffraction-less property of the Bessel-beam 
makes it a natural choice for imaging scattering specimens which are 
otherwise inaccessible by other means of microscopy imaging tech- 
niques such as, widefield, CLSM, 4PI, structured illumination micro- 
scopy and others. Notably, the experimental results show substantial 
minimization of the concentric ring effect, thereby facilitating a fine 
scanning ability without causing undue background fluorescence. 

Confocal laser scanning microscopy (CLSM), Structural illumina- 
tion microscopy (SIM) and other scanning-based fluorescence imaging 
techniques have advanced the fields of biophysics, biophotonics and 
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Figure 8 | Images of clusters of fluorescent polymer coated CaC03 particles (size 3-10 jum) inside a tissue-like gel matrix at depths of 108 /zm, and 
358 jum from the sample boundary. The location of this sample boundary is indicated by the red line and both positions 1 and 2 are located with respect to 
this line. A few out-of-focus clusters can be observed in position 1, one of these is encircled in blue. This same location is brought to focus by adjusting the 
detection arm along the optical axis (x-axis, see Figure 1.) and captured in position 2. A section in which a few coated CaC0 3 particles are clustured has 
been magnified. The scale bar (in white) represents 55 [im. 



nanoscale imaging. However, the inability to improve the penetration 
depth has limited the use of these techniques in thick scattering sam- 
ples such as tissue. The PSF of these imaging systems broadens at large 
depths in the sample, results in phase-mismatch and gives rise to poor- 
resolution. These effects rules out the very concept of high-resolution 
noise-free imaging. Widefield and CLSM imaging suffer from auto- 
fluorescence and off-layer (top and bottom of focal plane) excitation. 
With the improved depth penetration of Bessel-like beam, the illu- 
mination can be limited to just a pencil-like region of the specimen, 
thereby preventing premature photobleaching of the neighboring 
regions. Additionally, we observed a substantial reduction of the ring- 
ing effect surrounding the central pencil-like beam, which further 
reduce undue photobleaching and off-focus fluorescence. This illu- 
mination scheme coupled with orthogonal detection demonstrates 
the capability for selective imaging from a specific plane deep inside 
the specimen. In such a configuration, the lateral resolution of the 
illumination arm determines the axial resolution of the overall imaging 
system, which is beneficial for deep selective plane imaging. It is dif- 
ficult to go down to a resolution of few hundreds of nanometers 
because of the small working distance of high NA detection objective. 
Nevertheless, this configuration enables the study of thick scattering 
biological specimens at large penetration depths. 

The ability to look deeper inside a biological specimen has pro- 
found implications in biological systems. This ability is prohibited by 
the poor transmission of light through the specimen which com- 
prises several interfaces and scattering centers. To increase penetra- 
tion, one has to employ a beam that does not undergo significant 
diffraction. Specifically, the Gaussian beam is known to have signifi- 
cant diffraction. Bessel beams show much promise which is mainly 
because they are diffractionless and can self- reconstruct. When 
obstructed, the scattered field interferes with the unscattered field, 
resulting in constructive interference along the direction of beam 
propagation (which is same as the optical axis). Mathematically, 
the complex amplitude of the Bessel-like beam can be modeled as, 
U(r) = C 0 J 0 (kTp) e~i kzZ , for which the intensity distribution is, 



h = \U\ 2 = | Q | 2 ll(k T p), where, k T = ^Jk 2 x + k 2 y ; k x , k p & k z are the 
components of the wave-vector. The presence of a Bessel function in 
the intensity expression accounts for the side rings which decrease as 

p = \f(x 2 +y 2 ) increases and most importantly the intensity is inde- 
pendent of the propagation direction z, thus there is no spread of 
intensity. Therefore, the Bessel beam can propagate deep inside the 
scattering sample. In contrast, the intensity of the Gaussian beams is 



given by, I G = de 2 ? 2 ^ 2 ^; W(z) = W 0 



1 + 



1/2 



and has a in- 



dependence. Therefore the intensity spreads strongly at interfaces 
(refractive index mis-matches) in the biological specimen. 

The imaging technique has demonstrated quality imaging of both 
functionalized nanoparticles and polymercoated yeast in a thick gel 
matrix. PSF distortion is inevitable at large penetration depths espe- 
cially in a scattering sample. This technique demonstrates that Bessel 
beam-based excitation extends the imaging limit to almost 650 \im. In 
contrast, Gaussian beam-based excitation techniques are limited to 
~ 200 microns. Moreover, in an ideal isotropic medium, the computa- 
tional study revealed that the PSF retains its shape at greater penetra- 
tion depths. We have already witnessed in our gel-matrix sample that 
this does not hold true for z > 650 \im. In the computational study, 
scattering effects were not incorporated, which is a result of the dif- 
ficulty in modeling the scattering process in a complex biological 
environment. The PSF undergoes substantial distortion at greater pen- 
etration depths in a realistic sample, such as the tissue-like gel-matrix 
sample. Therefore, we conclude that reliable imaging is possible for z < 
650 \im. Within this range, the recorded images (see, figures 8 and 9) 
show the capability for high quality imaging. In the future, this tech- 
nique will attract applications in particle tracking deep inside tissue and 
optical injection in addition to fluorescence imaging applications. 

Methods 

Excitation and detection scheme. The flexibility for obtaining information from any 
desired layer of the specimen is of paramount importance in fluorescence imaging. 
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Figure 9 | Fluorescent polymer (PAH conjugated with NHS-Rhodamine dye) coated Yeast cells are uniformly distributed throughout the tissue-like 
gel sample. Several images were captured by translating the detection arm by 50 /urn; three of these images are shown here, each highlighting a prominent 
cluster of Yeast cells. These clusters were observed at depths of 296 (im (A), 478 fim (B), and 616 jim (C) from the coverslip. The inset shows a zoomed 
portion of a few clustered Yeast cells. (D), (E) and (F) are the intensity plots along the white line passing through the clusters in these inset figures of 
(A), (B) and (C). The FWHM values (obtained by averaging the Y and Z profile) are indicated in each intensity plot. It is evident that, the FWHM does not 
appreciably increase at greater penetration depths thereby supporting the fact that resolution deep inside the specimen is restored. The scale bar (in white) 
represents 25 /urn. 



Assuming linearly polarized light illumination and spatial filtering techniques, the 
resulting electric field at the geometrical focus is given by 41,42 , 



Kxc(u,v4) = \E(u,v,A6)\ 2 
= | J 0 1 2 | 2 cos 2 (« + \I 2 \ 2 +2 cos(20) Reloi; . 

where, the modified diffraction integrals over the aperture angle are, 



(1) 



I 0 (u,v) = J A0 M(*) (1 + cos 6) sin 6xJ Q 



h(u,v)= J M M(*) sin 2 ^ 



s^Oe^Ue, and 



~I 2 (u,v)= J A0 M(*) (l-cos0) sin 9} 2 



s^fle't^fW u, vand^c 



the longitudinal coordinate, transverse coordinates and azimuthal angle 
respectively 42 . The function, M(*) = H[9 — 9 X ] — H[9 — a] represents the spatial 
filter, and H[a — b] = {0, a ^ b; 1, a > b} is the Heaviside step function. 0 1 is the 
parameter controlling the transmission characteristics of the filter through the 
outermost annulus defined by A9 = a — 9 X (see figure 1). 

We employ confocal theta detection scheme, where the detection is performed in 
the orthogonal plane, as shown in figure 1. The transformation matrix representing 
x'l [0 0 -ll IV 
/ = 0 1 0 y . This detection 
z'\ \_l 0 0 J \_z_ 
scheme has the advantage of high resolution at long working distances. Since, the 
emission process is isotropic and the fluorescent light is randomly polarized, the 
detection PSF is given by, 



the orthogonal plane is given by, 



h det (x'y,z') = \E(x',y',z')\ 



lo (*',/,*') +2\h{x',y',z')\ + hix'JJ) 



(2) 



(3) 



where the integration J 0 ,i,2 = Jg" = 0 (■ ■ on me aperture-free detection objective is 
carried over the detection aperture angle a det . Overall, system PSF of the proposed 
imaging system is given by, 



h sys (x,y,z) = hui(x,y,z) x h det (x',y',z'). 



(4) 



For the computational study, the excitation wavelength was chosen as 532 nm, for 
which the fluorescence emission was collected at 560 nm. We have chosen to work at 
a sampling rate (well above the Nyquist sampling rate) of 30 nm along the entire x, y, 
& z axes. For the simulation, the number of pixels chosen along the lateral x, y axes 
were 128 X 128 and 128 of these planes were considered along z-axis. This gave a 
lateral (xy) and axial (z) range of about 3.84 X 3.84 fim 2 and 3.84 fim respectively. The 
integration over the modified diffraction integrals, J 0 > ^i and I 2 was performed over 
the transmission angle (A9). The detection PSF is similar to the excitation PSF except 
for the scaling introduced by the Stokes -shifted emission wavelength of 560 nm. The 
detection scheme consists of scanning along the z-axis to successively collect light 
from different layers of the 3D specimen. 

Generation and characterization of a bessel-like beam. To generate a Bessel beam, 
we used aperture engineering techniques. A special optical binary mask was employed 
that allows the light to pass through the outermost annular region (see, 
Supplementary 2). This structured light is then allowed to pass through the 
back- aperture of the objective lens. The spatial filter has a radius of 12 mm and is 
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Figure 10 | The gel-matrix containing yeast cells (coated with multiple layers of polymer [PAH] conjugated with fluorescent dye [NHS-Rhodamin] and 
PSS), is imaged by using a confocal microscopy. Several XY images were acquired to form a 3D stack. The farthest observable yeast cell cluster was 
observed at a distance of Az = 227 fim from the sample boundary, (a) XY-plane of the 3D stack at a depth of Az = 227 fim. The insert shows the zoomed 
portion of the Yeast cell cluster, (b) and (c) are the corresponding YZ and XZ planes, respectively, that are obtained from the 3D stack, (d) A schematic 
diagram depicting the location of the planes within the specimen and its distance from the coverslip. (e) For comparison, XY plane similar to (a) is also 
shown at a closer depth of 16 fim. (f) Magnified view of a small Yeast cluster taken from the highlighted portion of (e). (g) and (h) are the intensity plots 
along the white lines in the inset figures of (a) and (f), respectively, for comparison at two varying depths. The broadening of the FWHM at a penetration 
depth of 227 /urn is evident, indicating the deteriorating system resolution. Scale bars are inserted in the respective images. 



tailored in the workshop to fine tune the edges. The back- aperture of the objective lens 
(Olympus, UPlanFLN 10X/0.3NA) is approximately 15 mm and the light is 
transmitted through the outermost window of 12 - 15 mm. The characteristics of the 
spatial filter mask are detailed in Supplementary 2. Upon excitation and with the 
consequent focusing by the excitation objective, the light from the periphery 



undergoes constructive interference at the core along with low intensity rings due to 
further interferences. 

To characterize the beam, we placed a CCD camera (JenoptikMF coo /) directly at the 
geometrical focus of the illumination objective lens. The camera was then translated 
along the optical z— axis to collect the images of the XY-profile of the Bessel-like 
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beam. The translation was achieved using a micro -meter translator at an interval of 
10 fim within a range of 1 mm. The images were then stacked together to obtain the 
3D view of the field distribution as shown in figure 6. 

We employed theta detection in which the light was collected 90° to the optical axis. 
The detection objective (Olympus, UPlanFLN 10X/0.3NA or Olympus UMPlan FL 
20X/0.4NA) collected and collimated the light. This was followed by an iris that is 
introduced to reduce spherical aberration. The fluorescence light was then allowed to 
pass through a biconvex lens (focal length 150 mm) that focuses it on the CCD 
camera. The data were collection from different z-layers of the specimen by trans- 
lating the detection objective with incremental steps of 50 fim. To enable sharp 
focusing of the detection objective, the translator on which the sample was mounted 
was adjusted along the z-axis accordingly. 

The sample chamber was mounted on a holder that was placed on two ortho- 
gonally oriented micrometer translation stages (Holmarc ,TS - 50 - MwlO - 01, 
India). This arrangement allows the sample to be displaced with respect to the illu- 
minating beam both axially (z-axis) and along the transverse x-axis with a resolution 
of 10 fim. We used cover-glass to build a sample chamber which ensured that we have 
a perfect parallel plane facing the illumination beam and the detection objective. 
Having an irregular surface facing the detection objective increases spherical 
aberration and as a result produces image distortion. 

Sample preparation. Materials. Poly(styrenesulfonate sodium salt) (PSS; MW = 
70,000), poly(allylamine hydrochloride) (PAH; MW = 15,000), Ammonium 
bicarbonate(Mf 4 HC0 3 ), calcium nitrate(Ca(AT0 3 ) 2 ), Sodium chloride (NaCl) and 
dimethyl sulfoxide (DMSO) were obtained from Sigma- Aldrich. The fluorescent dye 
NHS-Rhodamine (MW: 528, Ex/Em wavelength: 552/575) and fluorescent dye 
removal columns were purchased from Thermo Scientific. Agar powder was acquired 
from Sd Fine Chemicals. An optical filter with a cutoff at 550 nm was purchased from 
Thorlabs for filtering the incident light. 

System PSF characterization. To characterize the Bessel beam, we prepared a fluor- 
escent gel sample. 100 mg of Agar powder was dissolved in 20 ml of distilled water 
and kept at 100°C for 10-15 min. The liquid gel was then allowed to cool down to 
40°C. Then, 500 filit of molten agar gel was added to 100 filit of NHS-Rhodamine dye 
with constant stirring. This mixure was then poured into a detachable cubical 
chamber and allowed to solidify. Special care was taken to ensure that smooth surfaces 
faced the excitation light as well as the detection camera to reduce undue lensing 
effect. Once solidified, the fluorescent -gel sample was ready for the PSF character- 
ization of the imaging system. 

Preparation of the PAH-NHS rhodamine conjugate. The labeling procedure 43 ' 44 con- 
sisted of mixing a PAH solution (1 mglml) prepared in double-distilled water, with 
the dye dissolved in DMSO (10 mM), in a 1 : 15 (w : w) ratio. The mixture was then 
placed in the incubator. Darkness was maintained overnight with continuous gentle 
stirring to avoiding photobleaching. Subsequently, the resulting PAH-NHS 
Rhodamine was purified using fluorescent dye removal columns. 

Coated yeast cell and CaC0 3 suspended gel samples. To prepare the gel sample, 100 mgoi 
Agar powder was dissolved in 20 ml distilled water and heated at 100°C for 10 
minutes. The melted agar gel was then cooled to 40° C. 500 filit of melted agar gel is 
added to 100 filit of coated CaC0 3 crystals with intense stirring. This was then shaped 
as a cube. The preparation of coated CaC0 3 is described in Supplementary 1. Special 
care was taken to obtain a smooth surface that included the surface facing the 
excitation Bessel-like beam and the adjacent side facing the detection objective. 

A similar technique was employed for preparing the Cell- suspended gel sample. 
Care was taken to cool the melted agar gel to 37°C. Then, 500filit of the gel solution 
was added to the coated yeast cells and mixed well by pipetting. The resulting mixture 
was then allowed to solidify in a cubical detachable container. 
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